THE SYMMETRIC AND UNIMODAL EXPANSION OF EULERIAN 
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Abstract. This paper was motivated by a conjecture of Branden (European J. Combin. 29 
(2008), no. 2, 514-531) about the divisibility of the coefficients in an expansion of generalized 
Eulerian polynomials, which implies the symmetric and unimodal property of the Eulerian 
numbers. We show that such a formula with the conjectured property can be derived from 
^i^S ' the combinatorial theory of continued fractions. We also discuss an analogous expansion 

for the corresponding formula for derangements and prove a (p, q)-analogue of the fact that 
the (-l)-evaluation of the enumerator polynomials of permutations (resp. derangements) by 
the number of excedances gives rise to tangent numbers (resp. secant numbers). The (p, g)- 
analogue unifies and generalizes our recent results (European J. Combin. 31 (2010), no. 7, 
1689-1705.) and that of Josuat- Verges (European J. Combin. 31 (2010), no. 7, 1892-1906). 
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O I 1- Introduction 

^^ , The Eulerian polynomials An[t) can be defined by 

I: EMt)^-, 



X 



°° x" l-t 



n=0 

Let &n be the set of permutations on [n] = {1, . . . , n}. For a = (t(1) . . . o"(n) € ©„, the entry 
C^ ' i G [n] is called an excedance (position) of o" if i < a{i) and the number of excedances of cr is 

J^ . denoted by exco". It is a classical result (cf. [FS70]) that the Eulerian polynomials have the 

QQ I following combinatorial interpretation 

in 



An{t) = Y, *''''^'' = Mn, 0) + A{n, l)t + ■ ■ ■ + A{n, n - I)*"" 



o-G6,i 



^^ ' with A(n, k) = #{(7 € ©n '■ exco" = k} and the expansion 

o ■ 



L(n-1)/2J 

An{t)= Yl anM'a + tr'''^^ (1) 

A:=0 

where a^^k ars nonnegative integers with known combinatorial interpretations. Recall that a 
C^ ' sequence of real numbers oq, ai, . . . , a^^ is said to be symmetric if Oj = a^-i for i = 0, . . . [d/2\ 

and siad to be unimodal if there exists an index 1 < j < d such that ao < oi < • • • < 
cbj-i < ttj > flj+i > • • • > Orf. Note that the expansion ([TJ enables to derive immedi- 
ately the symmetry and unimodality of the Eulerain numbers {A{n,k)}Q<k<n-i- Foata and 
Strehl |FS76] studied the combinatorial aspect of the expansion ([1]) via a group acting on 
the symmetric groups. In 2008, generalizing Foata and Strehl's action, Branden [B ra08] 
gave a (p, g)-refinement of ([T]) and conjectured that the corresponding polynomial coefficient 
a„^fc(p, g) has a factor (p + q)^ for all < A; < [(n — l)/2j (see [BraOBl Conjecture 10.3]). In 
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this paper we shall give a new approach to such an expansion with a proof of his conjecture 
as a bonus (cf . Theorem [2|) . 

Next we shall study the derangement counterpart of ([1]). Recall that a permutation a G ©„ 
is a derangement if it has no fixed points, i.e., a{i) ^ i for all i G [n]. Let 2)^ be the 
set of derangements in ©„. The derangement analogue of the Eulerian polynomials (see 
[F570l[B?i90] ) is defined by 

Bn{t) = Y^ t'''""' = B{n,l)t + B{n,2)t'^ + ■ ■ ■ + B{n,n - 1)^"-^ 

with B{n,k) = #{a G S)„ : exccr = k}. The generating function for Bn{x) reads as follows 
(see lEETOlEiienilKM]): 

/ J^n[t) r — —fz — -. 

^ — ' n! e — te 

n=0 

By analytical method, one can show (cf. |Bra061 fZhaOS] ) that there are non negative integers 
bn,k G ^ such that 

Ln/2J 

Bn{t)= Ybn^H^+tr-'"- (2) 

k=l 

However, no combinatorial interpretation for bn,k seems to be known hitherto. A special case 
of our results (cf. Corollary [9] and Theorem [TT]) will provide a combinatorial interpretation 
for the coefficients bn^k- 

Another interesting feature of ([T|) and ([2]) is that they generalize the well-known relation 
between the (-l)-evaluation of A„(t) and Bn{t) and the Euler numbers E^ (see [FS70J ): 

A2n{-1) =0, A2n+l{-l)= {-lTE2n+l, .oN 

B2n{-1) =(-l)"^2n, 52„+l(-l)= 0, ^-"^ 

where the Euler numbers {-E'n}n>o are defined by X]n>o -^"^ ~ tanx + secx. It follows that 
0'2n+i,ri = E2n+i and 62n,n = -£'2n- We will provc a (p, g)-analogue of the formulae ([3]) (cf. 
Theorem [T0|) , which unifies the recent results in J JVIOI ISZlOj . 

Our main tool is the combinatorial theory of continued fractions due to Flajolet |Fla80] 
and bijections due to Frangon-Viennot, Foata-Zeilberger and Biane, see |FV791 IFZQOt IFlaSOj 
IBia931 ICSZ971 ISZlOj between permutations and Motzkin paths. 

2. Definitions and main results 

Given a permutation a G ©„, the entry i G [n] is called a descent oi a \i i < n and 
a{i) > a{i + 1). Also the entry i G [n] is called a weak excedance (resp. drop) oi a Hi < a{i) 
(resp. i > cr{i)). Denote the number of descents (resp. weak excedances, drops) in a by 
deso" (resp. wexa, dropo"). It is well known |FS70j that the statistics des and exc have the 
same distribution on &n- Since drop a = exccr"^ for any a G &n, we see that the Eulerian 
polynomial has the following two interpretations 

An{t) = Y^ f^^'" = Y *'^'°'"'- (4) 
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The two above interpretations of Eulerian polynomials give rise then two possible exten- 
sions: one uses linear statistics and the other one uses cyclic statistics. We need some more 
definitions. 

Definition 1. For a € &n, let (t(0) = a{n + 1) = 0. Then any entry a{i) ( ? E [n]) can be 
classified according to one of the four cases: 

• a peak if a{i — 1) < (T{i) and a{i) > a{i + 1); 

• a valley if a{i — 1) > a(i) and a{i) < a{i + 1); 

• a double ascent if a{i — 1) < (j{i) and a{i) < a{i + 1); 

• a double descent if a{i — 1) > a{i) and a{i) > a{i + 1). 

Let peak* a (resp. valley* a, da* a, dd* a) denote the number of peaks (resp. valleys, double 
ascents, double descents) in a. Clearly we have peak* a = valley* a + 1. 

For a G (3„, the statistic (31-2) a (resp. (13-2) a ) is the number of pairs {i,j) such that 
'i < i < 3 1^ n and (T{i — 1) > a{j) > a{i) (resp. a{i — 1) < a{j) < cr{i)). Similarly, the 
statistic (2-13) a (resp. (2-31) a) is the number of pairs (i, j) such that 1 < i < j < n—1 and 
a{j -f 1) > a{i) > a{j) (resp. a{j + 1) < a{i) < cr{j)). Introduce the generalized Eulerian 
polynomial defined by 

An{p, q, t, U, V, W) := Y, P^'^'^^^ "9^^^"^^ -^desa^da* a^dd* a^valley* a_ (5) 

Let &n.k be the subset of permutations a S &n with exactly k valleys and without double 
descents. Define the polynomial 

an,,{p,q)= Y. P^'-''^^g^''-'^''. (6) 

Note that the coefficient a2n+i,nip, q) is the {p, g)-analogue of tangent number, i.e., E2n+iip, q) 
in [SZTOl eq. (14)]. 

Theorem 2. We have the expansion formula 

L(n-1)/2J 

An{p,q,t,u,v,w) = Y an,k{p,q){twf{u + vtT-^-'^^. (7) 

fc=0 
Moreover, for all < k < [(n — l)/2j , the following divisibility holds 

{p + qf \an,k{p.q)- (8) 

Remark. Branden [BraOSj used the convention (t(0) = a{n + 1) = n + 1 for the definition 
oi peak (resp. valleys, double ascent, double descent). Let peakg, valley g, das, ddg be the 
corresponding statistics. To see how the above theorem implies his conjecture we define the 
reverse and complement transformations on the permutations a E 6„ by 

(T I— )■ o"'" = a{n)a{n — 1) . . . (7(1), 

o- h^ cj'= = (n + 1 - a{l)){n + 1 - a{2)) . . . (n + 1 - (T{n)). 

Clearly the reverse-complement transformation a 1— > a''^ := (cr*^)^ satisfies 

(des, peak*, valley*, da*, dd*, (2-13), (31-2) )(t 

= (des, valley^, peaks, das, ddfi, (13-2), (2-31))o-''^ (9) 
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Thus, setting t = u = v = lm Theorem [21 we have 

L(n-1)/2J 
^ ^(13-2)^^(2-31) a^peaks a ^ ^ a„,fc(p, g)2"-l-2'=it;'=. (10) 

Comparing two coefficients of w^ in both sides of (|lUp . we recover Branden's result [BraOSl 
(5.1)]: 

a„,.(p,g)=2-"+l+2. ^ ^(13-2).^(2-31)., (11) 

peak^ TT^k 

Moreover ([5]) conffi'ms his conjecture |Bra081 Conjecture 10.3]. 

Definition 3. For a € ©„, a value x = a(i) {i G [n]) is called 

• a cyclic peak if i = (T^-'^(x) < x and x > (t(x); 

• a cyclic valley lii = cj~'^{x) > x and x < (t(x); 

• a double excedance ii i = a~^{x) < x and a; < a{x); 

• a double drop if i = a~^{x) > x and x > cr{x); 

• a /ixed pomi if x = (t(x). 

Let cpeako" (resp. cvalley cr, cda cr, cddu, fix a) be the number of cyclic peaks (resp. valleys, 
double excedances, double drops, fixed points) in a. 

For a permutation a G 6„ the crossing and nesting numbers are defined by 

croscT = # {{i,j) G [n] x [n] : {i < j < a{i) < a{j)) \J {i > j > a{i) > a{j))} , (12) 

nest a = # {{i,j) G [n] x [n] : {i < j < a{j) < a{i)) \J {i > j > a{j) > a{i))} . (13) 

As in |Cor07 ] . we can illustrate these statistics by a permutation diagram. Given a G S„, we 
put the numbers from 1 to n on a line and draw an edge from i to (T{i) above the line if i is 
a weak excedance, and below the line otherwise. For example, the permutation diagram of 
a = 3762154 = (3 7 g 2 i 5 I) ^^ as follows: 




a = 3762154 

Clearly, there are three crossings (1 < 2 < a{l) < a{2)), (1 < 3 < cr(l) < cr(3)), (7 > 
4 > o-(7) > cr(4)) and three nestings (2 < 3 < o-(3) < a{2)), (5 > 4 > cr(4) > o-(5)), 
(7 > 6 > cr(6) > cr(7)), thus croso" = nesto" = 3. 

Definition 4. Given a permutation a G ©„, let (t(0) = and cr(n +1) = n + 1. The 
corresponding number of peaks, valleys, double ascents, and double descents of permutation 
a G &n is denoted by peakcr, valley a, da a, and ddu, respectively. Moreover, a double ascent 
a{i) of cr (i G [n]) is said to be a foremaximum if o"(i) is a left-to-right maximum of o", i.e., 
(t(j) < (t(z) for all 1 < j < i. Denote the number of foremaxima of a by fmaxa. 
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For instance, da(42157368) = 3, but da* (42157368) = 2 and fmax(42157368) = 2. Note 
that by the above definition we have peak a = valley a for any a € ©„. 

Theorem 5. There is a bijection $ on (3„ such that for all a G S„ we have 

(nest, cros, drop, cda, cdd, cvalley, fix)(T 

= (2-31, 31-2, des, da — fmax, dd, vaney,fmax)<I>((T). 

Since wex = n — drop, we derive immediately from Theorems [2] and [5] the following {p, q)- 
analogue of ([1]). 

Corollary 6. We have 

L(n-1)/2J 

o-e6„ fc=o 

Consider the common enumerator polynomial 

Bn{p,q,t,U,V,W,y) : = ^ pncsta^crosa^drop<x^cdaa^cdd<ry^cvalleyayfix<T (^^5) 

= y^ p(2-31)o- (31-2) crudes cr^dacr-fmaxcr^ddo-^valleyo-fmaxcr /-j^g-, 

o-e6„ 
When u = and t = v = 1 we can write 

Bnip,q,l,0,l,w,y) := ^ bn,k,j{p,q)w''y'^, (17) 

where bn^kjiP: q) is a polynomial in p and q with non negative integral coefficients. In order 
to give a combinatorial interpretation for bn^kjiP^Q) we introduce some more definitions. Let 
^n,k,j denote the subset of all the permutations a £ &n with exactly k cyclic valleys, d fixed 
points, and without double excedance, and let 6* ■ ■ denote the subset of all permutations 
a € &n with exactly k valleys and j double ascents, which are all foremaxima. We derive 
from Theorem [5] the following result. 

Corollary 7. We have 

bnAd{p,q)= Y. p'^'^'^'q''"'''' = J2 p('-^^)'^<Z^3'"'^^ (18) 

In particular, when j = 0, we obtain 

K,kfi{p,q)= Y^ pnesta^cTosa ^ ^ ^(2-31) a^(31-2) <7_ (^^g) 

Recall (see [SZlOj) that a coderangement is a permutation without foremaximum. Let S)* 
be the subset of ©„ consisting of coderangements, that is, £•* = {a £ (3„ : fmaxo" = 0} . For 
example, we have D| = {2143,3142,3241,4123,4132,4213,4231,4312,4321}. Thus, D„,fc,o 
is the subset of derangements a € !!)« with exactly k cyclic valleys, and without double ex- 
cedance, and £>* fc is ^^^ subset of coderangements a € £>* with exactly k valleys and without 
double ascents. The following is our main result about the polynomial Bn{p,q,t,u,v,w,y). 
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Theorem 8. We have the following decomposition 

l{n-j)/2\ 

Bnip,q,t,u,v,w,y) = ^y^ ^ 6„,fcj(p, g)(H'=(gu + tt;)"-^-^^ (20) 

j=0 k=0 

For any a € ©n, since exco" = n — dropo", setting y = in Theorems [8] we obtain immedi- 
ately the fohowing {p, g)-analogue of ([2]) . 

Corollary 9. We have 

ln/2] 
Y^ pnesta^cvosa^eycca ^ ^ 5„,fc,o(P, Q)t'' {I + gt)""^^ (21) 

ctGDu fe=0 

Note that the coefficient b2n,nflip, q) is the {p, (7)-analogue of secant number, namely, E2n{p, q) 
in [SZlOl eq. (15)]. A permutation a € S„ is said to be a (falling) alternating permutation 
if (t(1) > o"(2), cr(2) < o"(3), (t(3) > o"(4), etc. Let 2ln be the set of (falling) alternating 
permutations on [n]. It is a folklore result that the cardinality of 2l„ is the Euler number En- 
We derive from Corollaries [6l [7] and [9] the following (p, g)-analogue of ([3]) . 

Theorem 10. For n > 1, we have 

if n is even, 

^i-l)^ E.Ga„P^^"^^^"5^^^"^^" «/™ ^s odd; 



\^ /j^-jwexo- nestcr c 



\wexa-„nest CT„croscr 

0-G6,, 



and 



—^ I if n is odd. 

Remark. The p = 1 case of Theorem 1101 was proved by Josuat- Verges [JVlOl Eq. (5)-(6)]. 
Since the inversion number of any a € &n can be written as drop a + cros a + 2 nest a (cf. 
|SZ10j ). the p = q^ case of Theorem [TU] was used by Shin and Zeng [SZlOj Eq. (12)-(13)] to 
derive a g-analogue of ([3]) using the inversion number g-analogue of Eulerian polynomials. 

For any permutation a G S„, we denote by eye a the number of its cycles. Define 

Cn,m= Yl /^'"''^ (22) 

o-gS)„(fc) 

where Sn(fc) is the subset of derangements in Xi^ with exactly k cyclic valleys and without 
cyclic double descents. Clearly CnflW) = for all n > 1. The following result gives another 
generalization of the expansion ([2]) . 

Theorem 11. We have 

Ln/2J 
Y^ ^cyca^exca ^ ^ c„,fc(/3)t^(l + t)"-^^ (23) 

CTGD,! fc = 

The rest of this paper is organized as follows. In Sections [3H6] we shall prove Theorem O 
Theorem [5l Theorem [8] and Theorem II 11 respectively. In Section [7| we give two variations of 
the unimodal and symmetric expansion of Eulerian polynomials and one for its derangement 
analogue. Finally we conclude with an open problem related to the descent polynomial of 
involutions. 
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3. Proof of Theorem [2] 

A Motzkin path of length n is a sequence of points (sq, • • • , s„) in the plan N x N such that 
So = (0, 0), Si — Si-i = (1, 0), (1, ±1) and s„ = (n, 0). Denote by Tin the set of Motzkin paths 
of length n > 1. We shall call a step (sj_i,Si) East, North-East, South-East, respectively, 
if Si — Si-i = (1,0), Si — Si-i = (1, 1), Si — Sj_i = (1, —1). The height of the step (sj_i, Si) 
is the ordinate of Sj_i. Given a Motzkin path 7, if we weight each East (resp. North-East, 
South-East) step of height i by Oj (resp. bi and q) and define the weight of 7 by the product 
of its step weights, denoted by 'w{'y), then 

1 + E E ^(^)^'' = ^ -2 — • (24) 

„>i 7eOT„ I -box ^^^i^ ^ 

1 — OiX 



It is convenient to use two kinds of horizontal steps, say, blue and red. Let 9Jt^ be the set of 
colored Motzkin paths of length n. A Laguerre history of length n is a couple (7, (pi, . . . ,Pn)), 
where 7 is a colored Motzkin path of length n and {pi, . . . ,pn) is a sequence such that 
< Pi < v{si-i,Si), where v{si-i,Si) = /c if Sj_i = (z — l,k). Denote by Sjn the set of 
Laguerre histories of length n. 

Let a € S„, the refinements of three generalized patterns are defined by 

(31-2)^ a = #{i:i + l<j and a{i + 1) < cr(j) =k < aii)} , 
(2-31);. cr = # {i : j < i - 1 and a{i) < a{j) = k < a{i - 1)} , 
(2-13)^ a = # {i : j < i - 1 and a{i - I) < a{j) = k< a{i)} . 

We clearly have (31-2) a = ELi (31-2)fcCT, (2-31) cj = ELi (2-31)^ (^, and (2-13) ct = 
Y12=i (2-13)fc<7. Two numbers 4 = (31-2)^0" and r^ = (2-31)^(7 are called the left embracing 
numbers and right embracing numbers of fc G [n] in cr. 

We need a modified version of Frangon-Viennot's bijection ^fv '■ Sn — ^ -^n-i defined 
by the following: for a G ©„ we construct the Laguerre history (sq, • • • , Sn-i,Pi, ■ ■ ■ ,Pn-i), 
where sq = (0,0) and the step (si_i,Sj) is North-East, South-East, East blue and East red 
if i is a valley, peak, double ascent, or double descent, respectively; while pi = (2-13)^ a for 
i = 1, . . . ,n— 1. If hi is the height of (si_i, Sj), i.e., Sj_i = {i—l, hi), then (2-13)^ o"+(31-2)^ a = 
hi. Since cr(0) = a{n + 1) = 0, so n must be a peak and valley a = peak cr — 1. Thus 
(so, • • • , Sn-i,Pi, ■ ■ ■ ,Pn-i) is a Lagucrre history of length n — 1 and 



W{a) = tER7+NE7^EB7^ER7^NE7 JJ , 



where NE7, EB7, and ER7 are the number of North-East steps, East blue steps, and East 
red steps of 7. Therefore, 



n-l 



An{p,q,t,u,v,w)= Yl t'^^^+^^^^^^^^^^^r^^^^nt^^ + lW (25) 

where [n]p^g = (p" — q^)/{p — q). Given a Motzkin path 7, weight each step at height k by 

ak ■■= tw[k + l\p^q, bk := {u + tv)[k + l]p^q, Cfc := [fc + l]p,q (26) 
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if the step is North-East, East, and South-East, respectively, and the weight of 7 is defined to 
be the product of the step weights. Then the last sum amounts to sum over all the Motzkin 
paths of length n — 1 with respect to ([26]) . i.e.. 



An{p,q,t,u,v,w) = ^ ti;(7). (27) 

7GOT„_i 

It follows from ()24p that the equation ()27p is equivalent to the following continued fraction 
expansion 

1 



y^^An{p,q,t,u,v,w)x 



n-l 



n>l I - (^u + tv)[l]p qX [l]p.g[2]p,gtWX 

l-iu + tv) [2]p,gx - [^]p-^t^]p-^^^^' 

(28) 
For < A; < [(n — l)/2j, let an^kiP: <l:t,u,v) be the coefficient oiw^ in An{p, q,t,u,v,w), i.e., 

An{p,q,t,u,v,w) = ^ an,k{p,q,t,u,v)w''. (29) 

fc=o 

Substituting x ^ (u+tv) ^^'^ ^ ^ w{u ^ ^^ ([28j) . we obtain 

>p v^ an,k{P,q,t,u,v) ^ „_;l ^ 1 



1 - [l]p,ga; 






1 _ rol ^ _ [3]p,q[4]p,gWX 

J- L-Jjp,?-^ 



(30) 
Since the right-hand side of the above identity is free of variables t, u, and v, the coefficient of 
w x"'~^ in the left-hand side is a polynomial in p and q with nonnegative integral coefficients. 
If we denote this coefficient by 

E, / X an,k{p,q,t,u,v) 

Pn,k[p,q) ■■- 



t'^iu + tvY- 



l-2fc' 



then Pn^k{p-,q) = an,A:(p, 9, 1, 1;0) = a„,A;(p, Q'jO, 1). On the other hand, comparing ([29|) and 
dSD, we see that an,k{p,q) = an,kip,q,0,l). Thus Pn,k{p-,q) = an,kip,q)- This proves Q. 
Finally, as (p + <?) | ['^Jp.giJT- + l]p,g for all n > 1, each w appears with a factor {p + g) in the 
right-hand side of ([30]) . and the polynomial Pn,k{p, q) is divisible by {p + q)^ . D 

4. Proof of Theorem [5] 

In our previous paper [SZlOl Theorem 8], a bijection $ : ©„ — >■ ©„ was constructed such 
that 

(2-31, 31-2, des, fmax)cr = (nest, cros, drop, fix)$(o") Vo" G S„. 
In this section we first recall the bijection $ and then show that the same bijection satisfies 

(da — fmax, dd, valley)(T = (cda, cdd, cvalley)$(a) Va € (3„. (31) 
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Let a = cr(l) . . . cr(n) be a permutation of [n], an inversion top number (resp. inversion 
bottom number) of a letter i in the word a is the number of occurrences of inversions of form 
(i,j) (resp {j,i)) in a. We now construct r = $((t) in such a way that 

(2-31);, a = nestfe r Vfc = 1, . . . , n. 

Given a permutation ct, we first construct two biwords, ( /,) and (5) , and then form the biword 

■* ^ ' by concatenating / and g, and /' and g', respectively. The word / is defined as 



J' a' 

the subword of descent bottoms in a, ordered increasingly, and g is defined as the subword 
of nondescent bottoms in a, also ordered increasingly. The word /' is the permutation on 
descent tops in a such that the inversion bottom number of each letter a in /' is the right 
embracing number of a in a. Similarly, the word g' is the permutation on nondescent tops in 
a such that the inversion top number of each letter b in g' is the right embracing number of b 
in a. Rearranging the columns of r', so that the bottom row is in increasing order, we obtain 
the permutation r = $(cr) as the top row of the rearranged bi-word. 

Example. Let (7 = 41279658 3, with right embracing numbers 1, 0, 0, 2, 0, 1, 1, 0, 0. Then 

7\ A3 56\ /g\^/2A7 
J') VS 4 6 9^ ' VW U 2 7 5 3^ 

I - ( f 5\_A3562478 9\ /249385716^ 
^ ~ V/' 5'/ U 4 6 9 1 2 7 5 3y ^ VI 2 3 4 5 6 7 8 9^ 
and thus $(cr) = T = 249385716. 

Let r = ^{a) G ©„. We enumerate the triple statistics (da a — fmax a, dd o", valley cr) . 

• If A; is a double ascent of a and not a foremaximum of o", then /c is a ascent top and 
also bottom. So k belongs to g and g' . By definition, (31-2)^0" > 0. Hence, the 
column (^) does not appear in (5), i.e., t^^(A;) < k < T{k). So /c is a cyclic double 

ascent of r. Conversely, if T^^{k) < k < T{k), then the column k appears in g and g' 
and (31-2)^(7 > 0. It implies that fc is a double ascent of a and not a foremaximum 
of a. 

• If A; is a double descent of a, then A; is a descent top and also bottom. So k belongs to 
/ and /', then T~^{k) > k > T{k). Hence A; is a cyclic double descent of r. Conversely, 
if T~^{k) > k > T{k), then the column A appears in / and /'. It implies that A is a 
double descent of a. 

• If A is a valley of a, then A is a descent bottom and also ascent top. So A belongs 
to / and g', then T~^{k) > A < t(A). Hence A is a cyclic valley of r. Conversely, if 
r~^(A) > A < t(A), then the column A appears in / and g' . It implies that A is a 
valley of a 

Thus (|31|) is established. The proof of Theorem [S] is then completed. D 

We illustrate $ : 64 — )• 64 with their statistics in Figure [TJ 

5. Proof of Theorem [8] 
Let a € Gn, the refinements of crossing and nesting are defined by 

crosfc (7 = # {i G [n] : (i < A < a{i) < a{k)) \J {i > k > a{i) > o-(A))} , 
nestfc (7 = # {i G [n] : (i < A < a{k) < a{i)) \/ {i > k > a{k) > a{i))} . 
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(7 e 64 


r = $(cr) 


dcscr 


(31-2) cr 


(2-31) cr 


da cr — fmax cr 


ddcr 


valley cr 


fmax cr 


dropr 


crosr 
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fixT 
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1234 
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2 
















4213 


3142 


2 


1 
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Figure 1. Illustration of $ on ©4 with their statistics 

We clearly have cros a = J2k=i crosfc a and nest a = Yl^=i nestfc a. 

Using Foata-Zeilberger's bijection^i?^ • &n —^ ^n, we construct the Laguerre history 
(sq, . . . ,Sn,Pi, ■ ■ ■ ,Pn), where sq = (0,0) and the step (sj_i,Sj) is North-East, South-East, 
East blue and East red if i is a cychc valley, cyclic peak, cyclic double ascent (or fixed point), 
or cyclic double descent, respectively; while pi = nestj a for i = 1, . . . ,n. Then, we have 



nestj a + crosj a 



hi, 
hi - 
hi, 
hi - 



if (si_i,Sj) 
if (sj„i,Sj) 
if (si_i,Si) 
if (sj_i,Si) 



is North-East; 
is South-East; 
is East blue; 
is East red. 



Thus (so, ■ ■ ■ ,Sn,Pi, ■ ■ ■ ,Pn) is a Laguerre history of length n and 

n 
){a) = tER7+NE7^EB7^ER7^NE7yEB*7^NE7-HEB7TTyj,^h,-l-pi^ 



W{ 



where NE7, EB7, and ER7 are the number of North-East steps. East blue steps, and East 
red steps of 7 and EB* 7 is the number of East blue steps whose height is equal to pi. Given 
a Motzkin path 7, weight each step at height h by 



Ofc := tw[h + 1 



p.-?' 



bk ■■= yp^ + {qu + tv)[h]p^q, Ck ■■ 



■ip,q 



(32) 
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if the step is North-East, East, and South-East, respectively, and the weight of 7 is defined to 
be the product of the step weights. Then the last sum amounts to sum over all the Motzkin 
paths of length n — 1 with respect to ([32]) . i.e.. 



Bn{p,q,t,u,v,w,y) = ^ '^(7). (33) 

7Gan„ 

By ()24p . we have the following continued fraction expansion: 

Y] Bnip, q, t, u, V, w, y)x" = 5 ' (^4) 

"To 1-60X "^^ ^ 

1 — OiX 

where at = tw[h + l]p,g, bh = yp'^ + {qu + tv)[h]p^q, and Ch = [h]p^q. 

For < A; < [(n — d)/2\ , let bn,k,d{p, Q, t, u, v) be the coefficient of w^y'^ in Bn{p, q, t, u, v, w, y), 
that is, 

n l(n-d)/2i 

Bnip,q,t,u,v,w,y) = '^ ^ bn,k,diP^'i,t,u,v)w^y'^. (35) 

d=0 fc=0 

Substituting X ^ (^^, u; ^ ^^^(^^, and y ^ (qu + tv)y in ® , we obtain 

„ V{n-d)/2\ 

y- v^ y- On,fc,d(p,g,n,^',t) ^.fc d^n 

n>Od=0 fc=0 ^^ ^ 

1 



1 - (y + [0]p,g)x 



[l]p_qu;j; 



l-{py + [l]p,q)x 



1 - (p2y + [2]p,g)x 






Since the right-hand side of the above identity is free of variables t, n, and v, the coefficient of 
w y x^ in the left-hand side is a polynomial in p and q with nonnegative integral coefficients. 
Denote this polynomial by Qn,k,d{p,q), then 

r, I \ bn,k,d{p,q,t,u,v) 

Qn,kAP,'i)--= t^^q^ + tvr-^k-d - (36) 

Hence Qn,k,d{p, q) = bn,k,dip, q, 1, 0, 1). By (^ and ([35]), we derive that Qn,k,d{p, q) = bn,k,d{p, q), 
which is given by (jlSp . This completes the proof. D 



Remark. Let An{p, q, t) = An{p, q, t, 1, 1, 1) and Bn{p, q, t) = Bn{p, q, t, 1, 1, 1, 1). From ([28]) 
and p^ we derive 

E Mp, 1, t)x- = 1 + ^'^^'^",1 .,1 , . , (37) 

1 - (1 + mU- r2] [31 tx^ 
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and 



^i3„(p,g,t)x" = ^^^^, . (38) 



2 +^2 



[2\„^tx 



1 - mp,, + [2].,,)^ - '-^^^^ 

In view of the following contraction formulae with C2i-i = [i]p,q and C2i = t[i]p,q for i > 1: 

1 ClX 

=1 + ^ 



CiX -, / N C2C32;2 

1 1 - (Cl + C2)X 



C2X -, / , X C4C5X2 
1 1 - (C3 + C4JX 



1 — CiX 



ClC2X^ 



. . C3C4X^ 
1 - (C2 + CsjX 



we see that 



Y, MP: Q: t)^" = E ^"(P' ^' *)^" = ' (39) 

n>0 n.>0 ]^ ^^^ 

C2X 



_ C3X 



where C2i = t[i]q and C2i-i = [i]q. Hence An{p,q,t) = Bn{p,q,t), i.e., the two triple sta- 
tistics (2-13, 31-2, des) and (2-31, 31-2, des) are equidistributed on 6„. Although a bijection 
showing the latter equidistribution can be constructed by combining the known bijections 
from permutations to Motzkin paths, a direct bijective proof of the above equidistribution is 
desired. 



6. Proof of Theorem [TT] 
Consider the polynomial 

Ln/2J 
Cni/3,t,U,V,w) := Y^ ^cyca^exca^cdaa^cdda^cvalloya ^ ^ Cn,kW ,t, U, v)wK (40) 

From [Zen93] we know that 



l + ^C„(/3,t,u,v,i/;)x" 



n>l 



1 - 0(tu + v)x ^— '-- ^ ^ 

1 - l{tU + V)X — '— ^ 2" 



-. (41) 



3(/3 + 2)twx^ 
1 - 2{tu + v)x - -^ 
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w{tu+v)'^ 



Substituting x ^ (tu+v) ^^'^ ^ ^ t — ™ (j^T|) , we get 

1 + EE y^^"rl ^^^" = 7^^-^-2 • (42) 

^^ ^-^ t^itu + vY-^^ l(/3 + 0)wx^ 

n>l fc=0 1 1^ 

2(/3 + l)wx2 
1 — X 



1-2. ='<^+,^)'"^; , 

J _ 3^. _ 4(/i + 3)«.i-' 



By the right-hand side we see that c* ^(/3) := iiTu^/Jn'^L must be a polynomial in /3, which 
is independent of u, v, and i. Hence we have 

Ln/2J 

a(/3,i,7i, .;,«;) = Y, clkm\tu + vr-^^w\ (43) 

fc=0 

Thus C„(/3, 1,1,0, w) = EK^ <,fc(/3)^^''- On the other side, by ([22]) and (gS]) we have 

C„(/3, 1,1,0,10) = Et=o^ c„,fc(/3)w;^ So c„,fe(/3) = <_;,(/3). Finally we get m by putting 
u = '(; = u; = lin (US]). ' D 



7. Two STAR VARIATIONS 

We show that the polynomial An, which is originally defined using linear statistics, has a 
new combinatorial interpretation by cyclic statistics. For cr = (t(1) . . . (j{n) G S„ and some 
integral function 'stat' on ©„, we define the star transformation from the permutation a to 
the function a* = cj*(l) . . . a*{n) from [n] to {0, . . . , n - 1} by 

a^a* = (ct(1) - 1) . . . (a(n) - 1). (44) 

For any statistic stat on a* we can define the corresponding star statistic 'stat*' on a by 
stat*((T) := stat ((7*). For instance, we can define the star cyclic statistics 

fix* o- = fixo-* = # {i G [n - 1] : i = a*{i) = a{i) - 1} , 
wex* a = wexiT* = # {i & [n — 1] : i < a*{i) = a{i) — 1} (= excu), 

cros* a = CTOsa* = # {{i,j) G [n] x [n] : {i < j < a*{i) < a*{j)) V {i > j > a*{i) > a*{j))} , 
nest* a = nesta* = # {{i,j) G [n] x [n] : {i < j < a*{j) < a*{i)) W {i > j > a*{j) > a*{i))} , 
cdd* fj = cdda* = # {i G [n - 1] : {(J*)-^{i) > i > a*{i)} , 
cda* a = cdacT* = # {i g [n - 1] : (cr*)"^(i) < i < a*{i)} , 
cvalley* a = cvalley ct* = # {i e [n - 1] : (cr*)"^(i) > i < a*{{)] . 

Clearly, a*{0) and {a*)^^{n) are neither defined nor needed. 

Remark. The star crossing number cros* (a) of a is different with the ordinary crossing 
number cros((T) of a. For example, for the previous example a = 3762154, we have a* = 
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-1 234567> 



2651043 =Q---)- 




a* = 2651043 

Since the four pairs (1 < 2 < a*(l) < a*(2)), (6 > 5 > cj*(6) > o-*(5)), (7 > 4 > a*{7) > 
(T*(4)), and (7 > 5 > (T*(7) > o"*(5)) are only crossings in a*, cros*(cr) = cros((T*) = 4. On 
the other hand, the three pairs (2 < 3 < a*{3) < a*{2)), (5 > 4 > a*(4) > a*(5)), and 
(7 > 6 > 0"*(6) > o"*(7)) are only nestings in a*, thus nest*((T) = nest(o"*) = 3. 

Theorem 12. The two hextuple statistics (2-13,31-2,des, da*, dd*, valley*) and (nest*,cros*, 
drop* —l,cda* + fix*, cdd*,cvalley*) are equidistrihuted on (3„. In other words, we have 

An{p, q, t, U, V,W)= Y^ p'^'^'^** '"9"°'^* '"t'^^P* <x-l^cda* a+fix* a^cdd* <x^cvalley* a_ (45^ 

Proof. Using the bijection <I> in Section H] (cf. Theorem 5), we can give a new bijection ^ on 
6n- Given a permutation o" = (t(1) . . . cr{n) € ©„, let o" = (<t(1) + 1) . . . (o"(n) + 1)(1) G ©n+i 
and consider the permutation r = $((t) E ©n+i- Since the last element (7(n + 1) of o" is 1, 
the first element r(l) of r should be n + 1. So r(2) . . . r(n + 1) is a permutation on [n] and 
we can define the permutation ^(o") := r(2) . . . T{n + 1) € ©„. By Theorem [SJ the bijection 
$ satisfies 

(2-31, 31-2, des)(T = (nest,cros,drop)<^((T). 
Since (2-13, 31-2, des)a = (2-31 + des -n, 31-2, des -1)(T and 

(nest*, cros*, drop* — 1)^((t) = (nest -|- drop — n, cros, drop —l)^{a), 

we have (2-13, 31-2, des)(T = (nest*, cros*, drop* — l)^(cr). Hence, assuming o"(0) = a{n + l) = 
0, the two hextuple statistics (2-13, 31-2, des, da*, dd*, valley*) and (nest*, cros*, drop* —1, cdd*, 
cda* -|- fix*, c valley*) are equidistrihuted on S^- □ 

Example. Given o" = 41279658 3, we define o" = 5238 10 76941 with right embracing 
numbers 1,1,1,2,0,1,1,0,0,0. Then 

7\ _ A 2 4 6 7 \ l'g\ _ /3 5 8 9 10^ 
J'J " U 9 5 7 loj ' [g'J " V2 3 8 6 1 

' - ( J" 9\ _ (I 2 4 6 7 3 5 8 9 10\ /lO 3 5 1 4 9 6 8 2 7 
^ ~ V/' 5'/ U 9 5 7 10 2 3 8 6 1 j ^ V 1 2 3 4 5 6 7 8 9 10^ 
and thus r = $(a) = 10 3 5 1 4 9 6 8 2 7. So, ^(ct) = r(2) . . . r(10) =35149682 7. 

Let ©* be the set of all sequences a* where a € S„, that is, S* := {a* : a G S„,} . By 
definition (jM]), for all o" € ©„, it holds that 

(fix*, nest*, cros*, wex*, cdd*, cda*, cvalley*)cr = (fix, nest, cros, wex, cdd, cda cvalley)cr*. (46) 

Note that, for a given a* G ©*, the entries and n are neither a cyclic valley nor a cyclic 
peak. We illustrate ^ : ©4 — > ©4 with their statistics in Figure [2j 
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r = *(cr) 


T* e&l 


descr 


(31-2) cr 


(2-13) a 
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Figure 2. Illustration of ^ on S4 with their statistics 



Given a a G ©„, the diagram of a* £ 6* consists of cycles i ^ cj*(i) — ^ • • • — ^ i with 
i G [n], and the path n — ?• a*{n) — t- • • • — ?• 0. Let eye* a be the number of cycles in the 
diagram of a* . For example, for the previous example a = 3762154 and a* = 2651043, we 
have cyc*(o") = cyc(o"*) = 2, since one cycle 1 ^ 2 — > 6 — ?• 4 and one path 7 ^ 3 — ?> 5 — ?• 
exist in a*. Let &nik) be the set of permutations a G ©„ with cvalley* a = k and cdd* a = 0. 
Introduce the polynomial 

dn,m= E /3^^^*'^-''^*^ (47) 

(Tee„(fc) 



Theorem 13. We have 

L(n-1)/2J 
Y^ ^cyc* <T-fix* a^cxca ^ ^ (i„,fe(/5)t'=(l + t)""^" 
o-e6„ k=0 

Moreover, for all k > 1, the polynomial dn^kil^) has a factor (/? -|- 1). 
Proof. Let 



■2k 



(48) 



Dn{(3,t,U,V,w) := Y^ ^cyc*a-fix*a^wex*a^cda*a+fix*a^ 



cdd* (T„,, cvalley* cr 



(49) 



o-eSn 
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By the same method in |Zen93| to comit the cycles, we obtain 
^Dn{l3,t,u,v,w)x''-'^ = 

n>l 

1 



1 - l(tu + v)x — '—— ^ ^ 

^ ^ , , 2(p + 2)twx^ 

1 - 2(tu + v)x '-— ^ ^ 

1 — 3[tu + v)x 

Define the polynomial dn,k{f!^,t,u,v) to be the coefficients of w'' in Dn{l3,t,u,v,w): 



(50) 



L(n-1)/2J 



Dni(3,t,U,V,w) = ^ dn,k{P,t, 



\ k 

u,v)w . 



fc=0 

(tu+v) 



Substituting x ^ (ty^+v) ^^'^ ^ ^ ^ " ^^ (f50|l , we get 



E E '"'^^^"Tl ^v-^ = ^^^ . (51) 

n>l fc=0 1 — X 



1 o^ 2(/3 + 2)t^x^ 

1 3x 3(/3 + 3)t/;x^ 



So d^i^{P) := ^k u^ly\rl-\-2k must be a polynomial in /3, since the right-hand side of (l5T]l is 
independent oi u, v and t. Hence we have 

L(n-1)/2J 

Dn{P,t,u,v,w)= Y. dlkif3)t\tu + vr-'-^'w'. (52) 

k=0 

To verify ()47p . taking t = u = 1 and f = in ()49p and ()52p we obtain 



L(n-1)/2J L{n-1)/2J 

D„(/3,l,l,0,«;)= ^ dn,kmw''= Yl <,kif^) 

fc=0 A:=0 



'w''. 



Thus dn,k{(3) = d*^kil^)-> ^"^^ ^^ S^^ (f^8]l by putting u = v = w = 1 in ([52]) . Clearly, for 
all A; > 1, the polynomial dn^k{P) has a factor (/3 + 1) because of the presence of the term 
(/3 + l)w at the second row of the continued fraction on the right side of (|5ip . D 

8. Concluding remark 
Consider the descent polynomial of involutions on [n]: 

/„(t) := Y t"^"" = Hn, 0) + I{n, l)i + • • • + /(n, n - l)t"-\ (53) 

where 3^ is the subset of involutions in (3„. The sequence {/(n, A;)}o<fc<n-i is known [GZ06] 
to be symmetric and unimodal. F. Brenti had conjectured that the sequence {/(re, fc)}o<fc<n-i 
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is log-concave, which was later disproved by Barnabei et al. |BFSd9] . Hence the best result 
one could expect for /^(i) is the following expansion 



L(n-1)/2J 

I„{t)= Yl an,kt\^ + tT~'~"" with a„,fcGN. 

k=0 



-2k 



(54) 



This is Conjecture 4.1 in [GZ06J . Unfortunately, the generating function of In{t) does not 
have a nice continued fraction expansion as that for An{t) or Bn{t). 
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The first values of an,k{p-, 



Appendix 
are given by anfi{p, q) = 1 for 1 < n < 5 and 



n\k 







1 2 


1 








2 








3 






p + q 


4 






{p + q){p + q + 2) 


5 




{P + 


q)[{p + q? + 2(p + g) + 3] {p + qfip" +pq + q^ + l) 



For J = 0, 1, 2 the first values of the polynomials bn,k,jip, q) are given as follows. 



where 



j=o 


k=0 


1 






2 3 


n= 1 













2 













3 













4 











{p + q? + l 


5 









(p 


+ q)^ + 2{p + q)^ + 2 


6 







(p+qy 


+ 2{p + qf + 3{p + q)^+3 56.3,o(p,g) 



j = l 


k=0 


1 


2 


n=l 


1 






2 









3 





p + 2 




4 





2p + 2 




5 





3p + 2 


b5.2,l{P^l) 



j=2 


k=0 


1 2 


n= 2 


1 




3 







4 





p^ + 2p + 3 


5 





3p2 + 4p + 3 


6 





6p2 + 6p + 3 66,2,2 (p,g) 



^6,3,o(p, q) = ip + qf + (1 - 2m)(p + q) + (2 + p'q'){p + qf + 1, 

fe5,2,i(p, q) = {p^ + 2p + 2)g2 + (2/ + 4p2 ^ 4p)^ ^ (p4 ^ 2p3 + 2p^ + 2p + 3), 
&6,2,2(p, q) = (/ + 2/ + 5p2 + 4p + 3)g2 + (2p5 + 4/ ^ ^q^s ^ 3^2 ^ g^^^ 
+ (p^ + 2p^ + 5/ + 4p^ + Qp^ + 6p + 6). 
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The first non-zero values of c„,fc(/?) (2 < n < 7) are given by the following table. 



n\k 


1 


2 


3 


2 


P 






3 


P 






4 


P 


/3(3/3 + 2) 




5 


P 


2/3(5/3 + 4) 




6 


P 


/3(25/3 + 22) 


/3(15/32 + 30/3 + 16) 


7 


P 


4/3(14/3 + 13) 


/3(105/32 + 238/3 + 136) 



The first values of dn,k{P') are given by dnfl{(3) = 1 for 1 < n < 7 and 



n\k 





1 


2 3 


1 
2 
3 
4 
5 
6 
7 




/3 + 1 
4(/3 + l) 
ll(/3 + l) 
26(/3 + l) 
57(/3 + l) 


(/3 + l)(3/3 + 5) 
(/3 + l)(25;3 + 43) 
10(/3 + 1)(13;3 + 23) ip + l)(15/32 + 60;3 + 61) 
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